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ABSTRACT 

We present a spectroscopic analysis of the metallicities, ages, and alpha- 
elements of the globular clusters (GCs) in the giant elliptical galaxy (gE) NGC 
4636 in the Virgo cluster. Line indices of the GCs are measured from the inte- 
grated spectra obtained with Faint Object Camera and Spectrograph (FOCAS) 
on the Subaru 8.2 m Telescope. We derive [Fe/H] values of 59 GCs based on the 
Brodie & Huchra method, and [Z/H], age, and [a/Fe] values of 33 GCs from the 
comparison of the Lick line indices with single stellar population models. The 
metallicity distribution of NGC 4636 GCs shows a hint of a bimodality with two 
peaks at [Fe/H] = — 1.23(cr = 0.32) and — 0.35(cr = 0.19). The age spread is 
large from 2 Gyr to 15 Gyr and the fraction of young GCs with age < 5 Gyr is 
about 27%. The [a/Fe] of the GCs shows a broad distribution with a mean value 
[a/Fe] ~ 0.14 dex. The dependence of these chemical properties on the galacto- 
centric radius is weak. We also derive the metallicities, ages, and [a/Fe] values 
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for the GCs in other nearby gEs (M87, M49, M60, NGC 5128, NGC 1399, and 
NGC 1407) from the line index data in the literature using the same methods as 
used for NGC 4636 GCs. The metallicity distribution of GCs in the combined 
sample of seven gEs including NGC 4636 is found to be bimodal, supported by 
the KMM test with a significance level of >99.9%. All these gEs harbor some 
young GCs with ages less than 5 Gyr. The mean age of the metal-rich GCs 
([Fe/H] > —0.9) is about 3 Gyr younger than that of the metal-poor GCs. The 
mean value of [a/Fe] of the gE GCs is smaller than that of the Milky Way GCs. 
We discuss these results in the context of GC formation in gEs. 

Subject headings: galaxies: abundances — galaxies: elliptical and lenticular, cD 
- galaxies: individual (NGC 4636, M87, M49, M60, NGC 5128, NGC 1399, 
NGC 1407) — galaxies: star clusters: general 



1. Introduction 



It is generally believed that the stars in a globular cluster (GC) are born almost at 
the same time and have similar metal abundances, ev en though there are some GCs with 



multiple stellar po pulations in the Milky Way (MW) flMarino et all 120091 : iLee et all 12009 



Bellini et al.ll2010l ). GCs are often assumed as single stellar population (SSP) systems and 
their observed quantities can be compared with theoretical models for the estimation of their 
physical parameters. They also contain information on the formation and evolution of their 
host galaxies. Therefore, analysis of their ages and metallicities can provide critical clues for 
understanding the evolution of their host galaxies. 

According to the current models of structure formation, giant elliptical galaxies (gEs) 
are considered to form and evolve through mergers with various galaxies. They might have 
accreted thousand s of GCs and sometimes formed new sta r clusters during the merging 
process (ILed 120031 ; iBrodie fc Straderl 120061 : ILee et al.l l2010al ). Therefore, investigation of 
ages and metallicities of the GCs in gEs may allow us to understand how GCs and their 
host galaxies formed and evolved. Though, in the case of gEs, it is not possible to obtain 
detailed information on the resolved stars in GCs with current telescopes, we can effectively 
estimate their age and metallicity using integrated spectroscopy. 

There are several studies on the ages and metallicities of GCs in nearby gEs based on 
integrated spectroscopy, which are summarized briefly below for six gEs (NGC 5128, M87, 
M 49, M60, NGC 14 07, NGC 1399). Basic properties of these galaxies are given in Table 3 
of ILee et al.l f|2010ah . 
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The GCs of the nearest g E NGC 5128 we r e extensively studied thanks to its rela- 



tively short distance, 3.8 Mpc flPeng et al 



(Peng et al. 


2004; 


Beaslev et al. 


2008: 



Woodley fc Gomez 2010 ). Recently, Woodley et al. ( 201dh derived the metallicities, ages 



and [a/Fe] values for 72 GCs in NGC 5128 using the Lick indices obtained with the Gem- 
ini/GMOS. Ninety-two percent of the metal-poor GCs (MP GCs) have ages older than 8 
Gyr, while 56 percent of the metal-rich GCs (MR GCs) do so. They found that the metal- 
licity distribution of the NGC 5128 GCs is bimodal if seen through the metallicities derived 
directly from the spectroscopic Lick indices ([MgFe]'), but the distribution is not well con- 
strained if they use Lick indices and SSP models. From these results together with the [a/Fe] 
data, they suggested that most GCs in NGC 5128 formed rapidly at early times, while some 
GCs formed later with subsequent major accretion. 

The GCs of M87 (NGC 4486), a c D galaxy located at the center of the Virgo cluster, 
were studied spectroscopically early by ICohen et al.l (119981 ). They reported that the mean 
metallicity of the GCs is higher than that of the MW GCs, and that the metallicity distribu- 
tion shows a bimodal feature with peaks at [Fe/H] = -1.3 and -0.7. The metallicities of the 
M87 GCs show a small radial gradient with a large scatter, but the ages of the GCs show 
no radial gradient. They suggested that the median age of the M87 GCs is 13 Gyr, similar 
to that of the MW GCs. 



M49 (NGC 4472) is the brightest gE in the Virgo cluster. ICohen et all (120031 ) presented 
the spectroscopic results of 47 GCs in this galaxy obtained with the Keck telescope. They 
concluded that the M49 GCs have slightly higher metallicities than the M87 GCs, while 
other metallicity and age parameters are basically similar to those of the M87 GCs. They 
estimated that the mean age of these GCs is about 10 Gyr. 



Pierce et al.l (120061 ) reported a study of 38 GCs in M60 (NGC 4649) in the Virgo cluster. 
They found that several young (2-3 Gyr old) GCs in this galaxy have supersolar metallic- 
ities, while the majority of the GCs have a wide range of metallicity ([Fe/H] = -2.0 to 0.0), 
and are older than 10 Gyr. They also found that the alpha-element ratios of the M60 GCs 
decrease with increasing metallicity. 

The GCs in NGC 1407, the brightest galaxy in th e NGC 1407 subgroup of the Eridanus 
group, were studied recently by ICenarro et al.l (120071 ). They found that 20 GCs are mostly 
old, while three GCs could be young (~ 4 Gyr). The metallicity of the GCs reaches up to 
slightly above solar and the mean [a/Fe] ratio is ~ 0.3 dex, similar to that of the MW GCs. 

Spectroscopic pro perties of the GCs in NGC 1399, gE in the cen t er of t he Fornax cluster, 
were studied early by iKissler-Patig et al.l (119981 ) and iForbes et al.l (120011 ) . From the Keck 
spectra of 10 GCs, IForbes et al.l (120011 ) showed that at least some of the GCs have supersolar 



-4- 



metal abundances. Two GCs with significantly higher H/3 values and enhanced abundance 
ratios could be regarded as young (~ 2 Gyr) or extremely old (~ 15 Gyr) ages with a 
warm blue horizontal branch. From this they proposed a possibility of a complicated age 
distribution among the MR GCs and of some constraints on their chemical enrichment at 
later epochs. 

These previous studies have some limitations. The number of the studied galaxies in 
this regard is only six, and the number of GCs in each galaxy is still small. The areal 
coverage for some galaxies is not wide enough, either. Therefore, it is needed to obtain 
new spectroscopic data for more GCs in each galaxy and for more gEs in order to better 
understand the chemical properties of gE GCs. 

We carried out a spectroscopic study of GCs in another Virgo gE NGC 4636. NGC 4636, 
located at th e southeast boun dary of the Virgo, harbors GC populations with a wide spatial 
distribution (ILee et al.ll2010bl ). It is one of the best targets for the study of age and metal 
abundance using the integrated spectroscopy together with M87, M49, and M60. There are 
several previous studies on the photometry, kinematics, and X-ray imaging of NGC 4636 
GCs (summarized below), but there is no published study of the metallicities, ages, and 
alpha-elements for the GCs based on spectroscopy. 

The specific frequency of NGC 463 6 GCs is SW = 5.6 — 8.9, larger than those of other 
elliptical galaxies in the Virgo Cluster (IKissler et al.lll994j ; iDirsch et al.ll2005l ). The GCs of 
NGC 4636 show a bimodal color distribution, while brighter GCs have intermediate colors 
(IDirsch et al.ll2005l ). The number of blue GCs is larger than that of red GCs. The number 
of the red GCs c onnected with low-mass X-ray binaries is ab out two times larger than that 
of the blue GCs kirn et al.ll20~0~6l ; IPosson-Brown et al.l[io09h . 

Using the radial velocity data for NGC 4636 GCs, ISchuberth et~al1 (bpOgh reported that 



the da rk matter fraction within one effective radius is 20% ~ 30%. However, IChakrabarty fc Raychaudhury 



(120081 ) sug gested tha t the d ark matter distribution is highly concentrated toward the in- 

) also confirmed the need of a large amount of dark matter. 



ner halo. 



Lee et al. (2010a 



Schuberth et al.l (120061 ) and lLee et al.l (j2010al ) found that the velocity dispersion for the blue 
GCs is slightly larger than that for the red GCs. They also studied the radial variation, the 
rotation, and the orbit of NGC 4636 GC sub-populations. 

We have been carrying out a project to investigate the spectroscopic properties of the 
GCs in nearby galaxies to understand the formation o f GCs in galaxies. Our study on the 



kinem atics of the GC system of M60 was presented in ILee et al.l (l2008bl ) and iHwang et al. 



(2008), an d another o n the GC system of the spiral galaxy M31 was given in 



Kim et al. 



(120071 ) and ILee et al.l (j2008d j. Recently, we presented the measurement of radial velocities 
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for the GCs in NGC 463 6 in a companion paper (IPark et al.ll2010l ). and detailed kinematic 
analysis of these data inlLee et al.l (j2010al ). We also presented a spectroscopic study of the 
M86 GCs ( TPark et alil2012bh . Here we present the chemical analysis of the NGC 4636 GCs. 
By combining the GC data of seven gEs (including NGC 4636) available in the literature, 
we also investigate the chemical properties of the gE GCs. 

This paper is composed as follows. Section 2 describes briefly the spectroscopic ob- 
servation and data reduction. In §3 we explain how we determine metallicities, ages, and 
alpha-elements of the GCs using the line indices measured in the optical spectra. §4 presents 
the chemical analysis results of NGC 4636 GCs. In §5, we compare the chemical properties of 
NGC 4636 GCs with those of other gEs, and discuss the implication of the results. Primary 
results are summarized in 86. 



Observation and Data Reduction 



We carried out a spectroscopic ob servation of GCs in NGC 4636, using the Faint Object 
Camera and Spectrograph (FOCAS; Kashikawa et al. 20021 ) at the Subaru 8.2 m Tele- 



scope. The spectra were obtained in Multi-Object Spectrography (MOS) mode of FOCAS. 
Target selection, spec troscopic observat ion, data reduction, and velocity measurement were 



described in detail in IPark et al.l ( 120101 ). We obtained flux-calibrated spectra of the targets 
using BD+33d2642, a standard star for flux calibration taken together with the targets. For 
the calibration of line indices, we used the spectra of five MW GCs (M5, M13, M92, M107, 
and NGC 6624) observed in the long slit mode during the same observing run. 

In Figure [TJ we show the flux-calibrated spectra of an MR GC and an MP GC in NGC 
4636, NGC 6624 and M13 in the MW Galaxy, and the NGC 4636 nucleus. These spectra are 
normalized at about 5400 A. These spectra show strong absorption lines typically seen in 
GCs such as H/3, Mgb, and Fel. MR GCs (Figure [T] (a) and (d)) have strong absorption lines 
at Mgb, while MP GCs (Figured] (b) and (e)) have relatively weak lines. Note that there 
are some sky-line residuals at ~ 5577 A, which are less than 2% of their original sky-lines. 
Although these residuals remain because of incomplete distortion correction in the optics, 
they affect little the index analysis in this study. 
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3. Metallicity, Age, and [a/Fe] Measurement 
3.1. Line index measurement 



We measured two 



rinds of absorption line indices from the spectra of NGC 4636 GCs: (1} 



the lBrodie fe Huchral (119901 ) line indices, used for determining the metallicities of GCs with 
empirical relations between absorption line indices and metallicities, and (2) the Lick line 
indices, used for measuring the metallicities, ages, and alpha-elements from the comparison 
with SSP models. Following subsections describe the measurement of each index. 



3.1.1. Brodie & Huchra (BH) line indices 



Brodie fe Huchral (119901 ) and iHuchra et al. I Jl996h presented linear relations between 



absorption line strength indices obtained from the integrated spectra of old stellar systems 
and their mean metallicity to determine the metal abundance of the systems. Their method 
was developed to minimize the systematic effects such as reddening, individual element 
abundance anomalies, and instrumental effects. They used the spectra and infrared colors 
of the MW GCs, M31 GCs, and NGC 188 giant stars for the metallicity calibration of the 
observed spectra. 

We measured the absorption line indices from the flux-calibrated spectr a of NGC 4636 
GCs a fter s hifting each spectrum to the rest frame following the prescription of lBrodie fe Huchra 
(119901 ) and lHuchra et al. I (119961 ) . The measured absorption line indices are calibrated to the 
BH index system with a zero point offset, Index(BR) = JncfexfSuba ru) + constant , deter - 
mined from the spectra of five MW GCs common in this study and IHuchra et al. I (119961 ) . 
The offsets (the constant in the above equation) we derived are 0.014 ± 0.015 for G band, 
-0.013 ± 0.009 for MgH, -0.021 ± 0.018 for Mg 2, and 0.008 ± 009 for Fe5270. Among 
the 105 NGC 4636 GCs with measured velocities (IPark et al.ll2010l ). we derived the BH line 
indices for 59 GCs with S/N > 10 and T\ < 21. The resulting index errors are smaller than 
0.10 mag for G band, 0.03 mag for MgH and Mg2, and 0.04 mag for Fe5270. 



3.1.2. Lick line indices 



Lick absorption line indices are useful for measuring the metall icity and age of old stellar 
systems from the integrated spectra with low/medium reso lution (IPuzia et al.l 120021 . l2005al ; 
Beasley et al.ll2008l ; iTrager et al.ll2008l ; IWoodley et al.ll2010l ). We mea sured the Lick line in- 
dices as follows. We first smoothed our spectra with the Lick resolution (jWorthey fe Ottaviani 
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19971 ) after shifting each spectrum to the rest frame. Then we derived the Lick line indices 



from the spect ra of NGC 4636 GCs, following t he definitions and the methods given in 



Wortheyl (119941 ) and IWorthey fc Ottavianil (119971 ). Line index errors are derived from the 



photon noise in the spectra before the flux calibration. Then we calibrated the resulting 
line indices to the Lick system with the zero point offset, Index (Lick) = ind ex (Subaru) + 



consta nt, d etermined from th e spect ra of five MW GCs common in this study, iTrager et al. 
f)l998h . and lKuntschner et al.l pjoj. 

Table [T] lists the zero point offsets derived in this study. The first, second, third, and 
fourth column are the name of the Lick index, unit of each index, zero point offset, and rms , 



respectively. Note that these offsets are slightly different from those in lPaudel et al. I ( 120101 ); 



these differences could systemically lead to the uncertainties of 2.3 Gyr in age, 0.14 dex in 
metallicity, and 0.18 dex in alpha-elements. We finally measure the Lick line indices for 33 
GCs in NGC 4636, listed in Tables |2] and [3j The associated errors are in Tables H] and [5j 
We use only the spectra with S/N > 15, where the S/N value is measured at the continuum 
part of the Lick line indices such as H/?, Mg 2 , Mgb, Fe5270, and Fe5335. 



3.2. Determination of Metallicity, Age, and [cn/Fe] 

We derived metallicity using three methods (the Brodie & Huchra method, the Lick 
index grid method and the x 2 minimization method), and age and [a/Fe] using two methods 
(the Lick index grid method and the x 2 minimization method) for NGC 4636 GCs. Each 
method is described below. 



3.2.1. Brodie & Huchra (BE) method 



We use the method described in iBrodie fc Huchra! (119901 ) to derive metallicities from 



the absorption line indices for NGC 4636 GCs. They recommended six best indices (G band, 
MgH, Mg2, Fe5270, CNB, and A) as the primary calibrators among the 12 line indices for 
the empirical relations between the BH line indices and the metallicities. We used, however, 
only four primary line indices to determine the metallicities of NGC 4636 GCs, excluding 
CNB and A because of low S/N values in their wavelength range. The relations between 
metallicity and the four indices are as follows. [Fe/H] Gkmd = 11.415 x Gband — 2.455, 
[Fe/H] Mffif = 20.578 x MgH - 1.840, [Fe/R] Mg2 = 9.921 x Mg2 - 2.212, and [Fe/H] Fe5270 = 
20.367 x Fe5270 — 2.086. We then derive the final metallicity value of each GC by taking 
an error-weighted average of four measurements. Here the error is the mean of the standard 
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deviation. The mean error for 59 GCs derived from this method is 0.31 ± 0.16 dex with the 
maximum error of about 0.6 dex. 



3.2.2. Lick index grid method 



To obtain the estimates of metallicity, age, and [a/Fe] from the integrated spectra of 
the NGC 4636 GCs, we use the Lick index grid method (the grid method hereafter), which is 



2005|). 



grids predicted from SSP models (ITripicco &: Belli Il995l ; iTrager et al. 



Traeer et al. 


2000l: 


1 


Thomas et al. 


iven by 


Thomas et al. 


( 


2003, 


2004, 



We use the grids of [MgFe]' versus H/3 provided by iThomas et al.l (120031 ) to estimate 
the metallicity and age of each GC. The composite index [MgFe]', defined as [MgFe]' = 
y/Mgb x (0.72 x Fe5 270 + 0.28 x Fe5335 ), is a good tracer of metallicity because of little 



sensitivity to [a/Fe] ( IThomas et al.l 120031 ). We use B.0 as an age indicator, which is the 
least sensitive to [a/Fe] among the Balmer lines (IThomas et al.ll2003l ). We also use the grids 
of Mg2 versus <Fe> in order to derive [a/Fe] estimates. The Mg2 index is very sensitive 
to [a/Fe], and the < Fe> index, defined a s <Fe> = (Fe5270+Fe5335)/2, is a metallicity 
indicator on this grid ( IThomas et al.l 120031 ). 



Figure [2] shows the measured indices in comparison with SSP model grids for various 
values of [Z/H] and age. Panel (a) displays the Lick line indices of H/3 versus [MgFe]' for the 



NGC 4636 GCs measured in this study as well as those for the NGC 5128 GCs (IWoodley et al. 
20101 ) for comparison. Grids indicate several values of [Z/H] (-2.25, -1.35, -0.33, 0.0, 0.35, 
and 0.67 dex) and ages (0.4, 0.6, 0.8, 1, 2, 3, 5, 8, 10, and 15 Gyr) from the SSP model 
with [a/Fe] = 0.2. The NGC 4636 GCs show a large dispersion in ages; most of them are 
old and some GCs are younger than 5 Gyr. The GCs in NGC 4636 and in NGC 5128 show 
similar age distributions, while the latter have a relatively small number of young GCs than 
the former. Metallicities of GCs in NGC 4636 and in NGC 5128 show a wide dispersion 
at [Z/H] < 0.7 dex. Panel (b) displays <Fe> versus Mg 2 , with the SSP grids for age = 8 
Gyr (mean age of our sample). While the scatter looks large, most GCs in NGC 4636 are 
concentrated on [a/Fe] of ~ and ~ 0.5. GCs in NGC 5128 also show a large dispersion 
and most of them are located at [a/Fe] ~ 0. In the followings, quantitative analyses on the 
age, metallicity, and [a/Fe] of NGC 4636 GCs are shown. 

Since the dependence of the H/3 versus [MgFe]' grid on [a/Fe] is not negligible, we need 
[a/Fe] information for accurate estimation of age and metallicity. An Mg 2 versus <Fe> grid 
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also needs an age value for accurate estimation of [a/Fe] and metallicity. Thus we need 
an iteration between the two grids in order to determine age, metallic ity, and [a/Fe] of a n 



object more accurately. We follow the iteration technique described in iPuzia et al.l (l2005al ). 
which is summarized in the following. First, the age and metallicity of one GC is determined 
from the H/3 versus [MgFe]' grid in a given [a/Fe] as in Figure [2] (a). Next, we determine 
the [a/Fe] and metallicity using the <Fe> versus Mg 2 grid with a previously determined 
age as in Figure [2] (b). Then, we derive again the age and metallicity of the GC using the 
H/3 versus [MgFe]' grid at the [a/Fe] value obtained above. We repeat this iteration until 
the [Z/H] difference converges to smaller than 0.01 dex. Finally, the converged values are 
selected as final values of age, metallicity, and [a/Fe] of the GC. 

To estimate the errors of the age, metallicity, and [a/Fe] values of a GC, we calculate 
ages, metallicities, and [a/Fe] values of the four data points composed of H/3 ± error and 
[MgFe]' ± error in the H/3 versus [MgFe]' grid and <Fe> ± error and Mg2 ± error in the 
<Fe> versus Mg2 grid. The difference between the average of these four values and the 
estimate directly calculated from the index is taken as the final error for the age, metallicity, 
and [a/Fe] of the GC. For some GCs outside the model grid, we adopt the value of the 
nearest enve l ope of the model grid in the direction of error vector as their estimates following 



Puzia et al.l (l2005al ). For example, if the data points are below the grid limit in the H/3 versus 
[MgFe]' diagram, we adopt the maximum value of age (i.e. 15 Gyr) at [Z/H] > —0.5 dex, 
but the age value of about 12 Gyr due to some overlapped grids at [Z/H] < —0.5 dex. For 
the data points below the grid limit in the <Fe> versus Mg 2 diagram, we also adopt the 
maximum value of [a/Fe] (i.e. 0.5 dex). 



3.2.3. x 2 minimization method 



Proctor et al.l (120041 ) suggested a method that can measure the physical parameters 
(metallicity, age, and alpha-elements) for GCs in a robust way based on simultaneous fitting 
of all the available Lick line indices with x 2 minimization. We call this method the x 2 
minimization method hereafter. This method maximizes the use of available information, 
while the grid method uses only a pair of indices. Recently this method was used in the 
st udies of the GCs in gE M60 by iPierce et al.l (120061 ) and early- type dwarf galaxies in Virgo 
bv lPaudel et~al~l (boioh . 



We applied this method to the NGC 4636 GCs using the same SSP models as used 
in the grid method. Among 10 Lick line indices except Mgi at 4500 — 5750 A, we used 
about eight indices after ~ 2a clipping of their x- The error of each parameter is the value 
where x 2 difference above the minimum reaches la significance level. We list the measured 
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metallicities, ages, and [a/Fe] for the NGC 4636 GCs in Table [6j 



4. Results 

4.1. A Catalog of Metallicity, Age, and [a/Fe] 



Among the 105 GCs whose radial velocities are measured in iPark et al.l (120101 ). we 
obtained [Z/H], age, and [a/Fe] for 33 GCs with the grid method and x 2 minimization 
method, while we obtained [Fe/H] for 59 GCs with the BH method. Table[6]lists metallicities, 
ages, and [a/Fe] values for the NGC 4636 GCs. 

Figure E] shows a comparison of the results for NGC 4636 GCs derived from three 
methods. Panel (a) shows a comparison of the metallicities obtained from the BH method 
and the grid method. [Fe/H] (grid or x 2 ) is converted from the total metallicity, [Z/H], us - 



ing the correlation between the two: [Fe/H] = [Z/H] — 0.94 [a/Fe] ( jThomas et al.ll2003l ). 
Here we adopt [a/Fe] = 0.2, which is the mean [a/Fe] for NGC 4636 GCs. Panel (a) 
shows a good agreement between the two measurements. We derive a linear least-squares 
fit: [Fe/H](BH) = 0.935(±0.080) [Fe/H] (grid) + 0.170(±0.069) with rms = 0.227. In the 
metal poor region, however, the metallicities derived from the grid method for three GCs 
are smaller than those from the BH method. Panel (b) illustrates a comparison of the 
metallicities obtained from the BH method and the \ 2 minimization method, showing 
a good agreement between the two measurements. We derive a linear least-squares fit: 
[Fe/H](BH) = 0.988(±0.072)[Fe/H](x 2 ) - 0.105(±0.078) with rms = 0.297. 

Panel (c) displays a relation between [Z/H] derived from the B./3 versus [MgFe]' grid 
and that from the <Fe> versus Mg 2 grid. It also shows a good agreement between the 
two estimates except for three outliers. We derive a linear least-squares fit: [Z/H](Mg 2 ) = 
0.914(±0.058)[Z/H]([MgFe]') +0.105(±0.031) with rms = 0.192. For the following analysis, 
we mainly use the metallicities of GCs measured from the H/3 versus [MgFe]' grid because 
the main concerns in this study are the ages rather than the alpha-elements of GCs. 

Panel (d) illustrates a comparison of the metallicities obtained from the grid method and 
the x 2 minimization method, which shows a good agreement between the two measurements. 
We derive a linear least-squares fit: [Z/H](x 2 ) = 0.946(±0.073) [Z/H] (grid) - 0.092(±0.052) 
with rms = 0.209. Panel (e) shows a comparison of the ages obtained from the grid method 
and the x 2 minimization method. Panel (e) shows a reasonable agreement between the 
two measurements with a large scatter. We derive a linear least-squares fit: Age(% 2 ) = 
1.017(±0.085)Age(grid) + 0.057(±0.739) with rms = 2.942. Note that there are six GCs 
with ages younger than 5 Gyr derived from both methods. The five out of six GCs have 
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relatively higher S/N (> 20). Therefore, the existence of young GCs is not simply due to 
the low S/N spectra (see next section for detailed discussion). Panel (f) shows a comparison 
of the [a/Fe] obtained from the grid method and the x 2 minimization method. It shows 
a reasonable agreement between the two measurements with a large scatter. We derive 
a linear least-squares fit: [a/Fe] (x 2 ) = 0.934(±0.064)[a/Fe] (grid) - 0.042(±0.028) with 
rms = 0.125. 



4.2. Distributions of Metallicity, Age and [a/Fe] 

Figure H] represents the distributions of metallicity, age, and [a/Fe] for the GCs in 
NGC 4636 derived from the three methods. Panel (a) shows that the [Fe/H] obtained from 
the BH method has a wide distribution at —2.2 < [Fe/H] < 0.4 with the mean value of 
[Fe/H] = -0.90 ± 0.12. Two peaks appear to be at [Fe/H] ~ -1.5 and - -0.4. The [Z/H] 
distribution obtained from the [MgFe]' grid method (Figure H](b)) shows a similar wide range 
of —2.4 < [Z/H] < 0.8 with the mean value of — 0.58±0.75. Two peaks are seen at —1.0 and 
—0.2. The results from the x 2 minimization method (dot-dashed line) are similar to those 
from the grid method (solid line). 

We performed the KMM test to analyze quantitatively the bimodality of the metallicity 



distrib ution of NGC 4636 GCs. The KMM test based on the algorithm of lAshman et al. 



(119941 ) gives an estimate of the improvement of the n-group fit over a single Gaussian. In our 
sample, the hypothesis of a unimodal metallicity distribution can be rejected at the 97.1%, 
71.4%, 98.4%, and 95.2% confidence level in the case of the [Fe/H] obtained from the BH 
method, the [MgFe]' grid, the color-metallicity relation, and the x 2 minimization method, 
respectively (see Figure E]). These suggest that the metallicity distribution of NGC 4636 
GCs is marginally bimodal. 

Figure [5] represents the metallicity distributions and their Gaussian fits determined 
from the KMM test. Panel (a) shows that two peaks are at [Fe/H] = — 1.23(<r = 0.32), 
and — 0.35(<7 = 0.19). The distribution of [Fe/H] determined with [MgFe]' grid (panel (b)) 
shows two peaks at [Fe/H] = -l.U(a = 0.22) and -0.41 (ct = 0.21). The distribution of 
[Fe/H] from the x 2 minimization method (panel (d)) also shows two peaks at [Fe/H] = 
— 1.22(cr = 0.25) and — 0.44(cr = 0.16). In addition we d erived the photometric metallicity 



distribut ion from the (C — Ti) color (IPark et al.ll2012al ) using the double linear equation 



given in iLee et al. I (j2008al ) (panel (c)). It shows two peaks at [Fe/H] = — 1.35(cr = 0.16) 
and — 0.59(cr = 0.28). We also plot the metallicity distribution of the MW GCs from lHarris 



(119961 ) (2010 version) in panel (a) for comparison. It also shows two peaks at [Fe/H] 



— 1.52(cr = 0.39) and — 0.52(cr = 0.23), each of which is at slightly lower values than the 
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NGC 4636 GCs. 

Figure H] (c) shows that the age distribution of the NGC 4636 GCs has a wide range 
from young (2 Gyr) to old (15 Gyr) with the mean grid age of 7.9 ± 0.8 Gyr and the mean 
X 2 fit age of 8.7 ± 2.0 Gyr. The results from both methods are similar. It seems that there 
are three peaks at 3, 8, and 14 Gyr. If we divide the age distribution into three bins of t < 5 
Gyr, 5 < t < 10 Gyr, and t > 10 Gyr, the number ratios are 9:13:11 (27:39:34%) in the grid 
age and 9:10:14 (27:30:43%) in the \ 2 fit age. The fraction of the young GCs with age < 5 
Gyr in our sample is about 27%, which is a slightly larger than that (18%) in NGC 5128 



( IWoodley et al 



20101 ) 



To check the robustness of the large fraction of young GCs in NGC 4636, we first divide 
the GCs into two groups based on their S/N: 22 GCs with S/N > 20 and 11 GCs with S/N 
< 20. The GCs in the two groups do not show any difference in their age distributions. 
There are six young (age < 5 Gyrs) GCs among 22 GCs with S/N > 20. Similarly, there 
are three young GCs among 11 GCs with S/N < 20. These suggest that the existence of 
young GCs is not simply because of spectra with low S/N. However, Brodie et al. (2005) 
argued that young or intermediate GCs may be due to relatively low S/N spectra or to sky 
subtraction difficulties. We also computed the probability of the existence of young GC with 
age < 5 Gyr and S/N > 20 using a simple simulation with GCs of age = 9 Gyr (a mean age 
of NGC 4636 GCs) and of Hbeta error = 0.34 A (a mean value of measured errors). The 
fraction of young GCs with age < 5 Gyr derived from simulations is 12%, which is smaller 
than the observational fraction (27%). This result indicates that the young GCs in NGC 
4636 could be caused in part by observational uncertainty. The reliability of the existence 
of young GCs needs to be checked with better data in the future. 

Alpha-element abundances indicate the duration of star formation in a GC. The GCs 
with an enhanced [a/Fe] ratio experienced rapid star formation over less than 1 Gyr timescale , 



while solar [a/Fe] GCs have an extended star formation history ( IBrodie fc Straderl 120061 ). 
Figure H] (d) displays the distribution of [a/Fe]. The overall distribution in the grid method 
shows a single peak with a large dispersion except for those at both ends, while that from 
the x 2 minimization method does not show a dominant peak. The mean [a/Fe] is derived 
to be 0.20 ± 0.11 in the grid method and 0.07 ± 0.05 in the x 2 minimization method. This 
mean value is between the values for the NGC 5128 GCs ([a/Fe] = 0.13 ± 0.03 in the grid 
method and [a/Fe] = 0.06 ± 0.03 in the x 2 minimization method, see Table ITU]) and for the 
MW GCs ([a/Fe] = 0.36 ±0.01). This suggests that the formation period of NGC 4636 GCs 
might have been more extended, as in the case of NGC 5128 GCs, than the MW GCs. 
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4.3. Radial distribution 

We investigated the distributions of the metallicity, age, and [a/Fe] of the NGC 4636 
GCs as a function of the galactocentric distance. Figure [6] (a), (b), and (e) show that [Fe/H] 
and [Z/H] (from BH, grid, and x 2 minimization methods) vary little with the galactocentric 
radii at R < 8' (< 34 kpc), but show a large scatter. It is noted that four GCs at 8' < 
R < 10'. 5 have on average ~ 0.5 dex lower [Fe/H] values than the GCs in the inner region. 
Figure [6] (c) and (f) show that the mean age of the NGC 4636 GCs changes little with the 
galac tocentric radii and the dispersion in age is large. This trend is similarly found in M60 



GCs (IPierce et al.ll2006l ). Panels (d) and (g) show no radial dependence of [a/Fe]. 



4.4. Relation between metallicity and age 



Figure [7J displays the relations between metallicity and age for the NGC 4636 GCs in 
comparison with NGC 5128 GCs and the MW GCs. The NGC 4636 GCs show, on average, 
an anti-correlation between age and metallicity. Younger GCs in NGC 4636 have higher 
metal abundances in bot h grid and y 2 minim ization results. It is also seen in the case of 
NGC 5128 GCs (see also lWoodlev et all (boioh and the MW GCs. This c orrelation is also 



seen for the GCs in M60 JPierce et al.l l2006h and M86 JPark et alil2012bh . Note that this 



correlation could b e caused by the correlated errors of age and metallicity 



see the error 



boxes in Figure [7]) (ITrager et al. Il2000t iKuntschner et al. 



biased sample selection (e.g., extremely bright GCs) (jWoodley et al. 



2001 



Paudel et al. 



2010( 1 or by a 



2010)- The age of the 



youngest GCs in the MW is 7 Gyr, but NGC 4636 and NGC 5128 have a significant number 
of GCs younger than 7 Gyr. The relative fraction of old (> 10 Gyr) GCs in NGC 4636 is 
smaller than that of NGC 5128 GCs and much smaller than that of MW GCs. 



5. Discussion 



In this study, we derived the metallicity, age, and [a/Fe] for the GCs in NGC 4636 
from the spectroscopic data. Here we compare these with those of GCs in other nearby 
gEs and discuss the results in the context of GC formation in gEs. We found six gEs for 
which the spectroscopic data of their GCs are available in the literature. We compiled 



the Li ck line index meas urements for the GCs in each gal axy from the 



ences: 



Cohen et al.l (119981) for M87, iBeaslev et al l fl2000h and lCohen et all (120031 ) for M 49 



Woodlev et all (l2010h for NGC 5128 



Pierce et al. 



fh006h for M60, 



ollow ing refer- 



Cenarro et al 



(120070 for 



NGC 1407, and iForbes et al.l (120011 ) for NGC 1399. To have a fair comparison of physical 
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parameters among gEs, we re-derived the metallicity, age, and [ct/Fe] for the GCs in each 
gE from the line indices using both grid and x 2 minimization methods as done for NGC 
4636 GCs (see Section 3.2). We also co mpared the G Cs in gEs with the MW GCs. The 
data of the MW GCs are obta ined from Harris! (1996 ) (2010 version ) for metallicity , from 
Marin- Franch et all (120091 ) and iDotter et al. (boioh for age, and from ICarretta et all (120101 ) 



for [a/Fe]. 



5.1. Metallicity distribution of GCs in gEs 



The color bimodality in recent ex tragalactic GC studies is one of the important in- 
gredients for GC formation scenarios (IBrodie fc Straderl 120061 ) . This bimodal color dis- 



tribution is interpreted as an intrinsic property due to a bimodal metallicity distri 


Dution 


(Brodie & Strader 


2006; 


Spitler et al. 


2008; 


Woodlev et al. 


2010; 


Chies-Santos et al. 


2012), 


or an apparent feature simply due to a nonlinear color-metallicity relation ( 


Yoon et al 


2006. 



2011alj bl). It is not yet clear which of these two is more supported by observational results. 
Here we investigate the metallicity distributions of the GCs in seven gEs. 

Figure [8] displays the distributions of metallicity derived with the grid method for the 
GCs in seven gEs. The metallicity distributions of GCs in most gEs seem to be bimodal as 
in the case of the MW GC system (see Figure (a)). We run the KMM test for quantitative 
analysis of the metallicity distribution, and list the results for the Gaussian peaks and their 
confidence levels in Table [71 There is no strong evidence for the bimodality for NGC 4636, 
M60, NGC 1407, and NGC 1399 that have a relatively small number of GCs. However, the 
bimodality for NGC 5128, M87, and M49 is supported by the KMM test with a confidence 
level > 98%. 

We combine the data derived from the grid method for all the GCs in seven gEs, and 
show the resulting metallicity distribution in Figure [8] (h). The combined gE sample (304 
GCs) shows clearly a bimodal distribution (155 MP GCs and 149 MR GCs). The KMM test 
yields two Gaussian peaks at [Fe/H] = -1.21(a = 0.32), and -0.42(<x = 0.27) with 99.9% 
confidence level. The MR peak for the combined gE sample is ~ 0.1 dex higher than that 
for the MW GCs ( [Fe/H] = -0.52(a = 0.23) ), while the MP peak is ~ 0.3 dex higher than 
that of the MW GCs ([Fe/H] = — 1.52(<r = 0.39)). It is noted that there are two more minor 
components at both ends ([Fe/H] < —2.0 and [Fe/H] > 0.0) in the combined gE sample. 

Figure [9] displays the distributions of metallicity derived with the \ 2 method for the 
GCs in seven gEs. The KMM test results for the Gaussian peaks and their confidence 
levels are given in Table [7J There is no strong evidence for the bimodality only for two 
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gEs (NGC 1399 and NGC 1407). However, the bimodality for NGC 4636, NGC 5128, 
M60, M87, and M49 is supported by the KMM test with a confidence level > 95%. The 
combined gE sample (315 GCs) shows clearly a bimodal distribution (167 MP GCs and 148 
MR GCs). The KMM test yields two Gaussian peaks at [Fe/H] = -1.25(a = 0.32), and 
— 0.42(cr = 0.25) with 99.9% confidence level, which are very similar to those from the grid 
method, [Fe/H] = -1.21(a = 0.32), and -0.42(cr = 0.27). It is noted that there are two 
more minor components at both ends ([Fe/H] < —2.0 and [Fe/H] > 0.0) in the combined gE 
sample, as found in the results from the grid method. 

To examine any age dependence of the metallicity distribution of the GCs in gEs, we 
divide the GC samples into two groups: young GCs with age < 10 Gyr and old GCs with age 
> 10 Gyr. In Figures |8] and EH we plot the metallicity distribution of the young GCs and old 
GCs (the dotted and dot-dashed lines, respectively). We also derive the mean metallicities 
of the young GCs and old GCs, which are marked with the vertical lines in Figures [8] and [9 



see a lso Table [8]). The mean values are determined with the biweight location (IBeers et al.l 



19901 ). and their errors indicate the 68% confidence levels obtained with a bootstrip method. 
The mean metallicity of all GCs in gEs ([Fe/H] = -0.80 ± 0.04 and -0.86 ± 0.04 from the 
grid method and the x 2 minimization method, respectively) is much higher than that for 
the MW GCs ([Fe/H] = -1.28 ± 0.05). The mean metallicity of the young GCs in each gE 
and the MW is 0.2 -1.3 dex higher than that of the old GCs. For the combined sample, 
the mean metallicities of the young GCs and old GCs are [Fe/H] = —0.55 ± 0.05 and [Fe/H] 
= -0.96 ± 0.05 in the case of the grid method, and [Fe/H] = -0.65 ± 0.06 and [Fe/H] 
= —1.00 ± 0.04 in the case of the \ 2 minimization method. Thus the metallicity of the 
young GCs is, on average, ~ 0.4 dex higher than that of the old GCs. It is noted that both 
old and young GCs appear to show a bimodal distribution. However, the MP component is 
dominant for the old GCs, while the MR component is dominant for the young GCs. 

We also examine a relation between metallicity and [ct/Fe]. We divide the sample into 
two groups: high [a/Fe] > 0.2 and low [a/Fe] < 0.2. The results are listed in Table |HJ The 
mean metallicity of GCs does not show any significant dependence on [a/Fe], except for 
two gEs (NGC 1407 and NGC 1399) with a small number of GCs. Note that some of the 
differences between galaxies could be dominated by selection effects (e.g., different brightness 
of GCs) or by small number statistics. 



5.2. Age distribution of GCs in gEs 



Figures [10] and [IT] show the age distributions of the GCs in gEs derived from the grid 
method and from the \ 2 minimization method, respectively. We also plot the age distribution 
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of the GCs in the combined sample for gEs (panel (h)). We divide the sample into subgroups 
according to their metallicity and [a/Fe]: MP groups with [Fe/H] < —0.9 and MR groups 
with [Fe/H] > —0.9, and low [a/Fe] groups with [a/Fe] < 0.2 and high [a/Fe] groups with 
[a/Fe] > 0.2. We derive the mean ages of each group, which are listed in Table |9J The mean 
values of the metallicity subgroups are shown by the vertical lines in Figures [10] and [TT1 

Several features are noted in these figures. First, all gEs in this sample show a wide 
range of ages for GCs, from 1 Gyr to 15 Gyr. This is in stark contrast to the case of the 
MW GCs. The MW GCs are all older than 7 Gyr and mostly older than 10 Gyrs (see Figure 
[7]). Second, the mean age for the MP GCs in the combined sample (12.2 ±0.1 Gyr from 
both methods) is about 3 Gyr larger than that for the MR GCs (9.2 ± 0.7 Gyr from the 
grid method and 9.7 ± 0.7 Gyr from the x 2 minimization method, respectively). The MW 
GCs show similar trends, but with a smaller difference of 1 Gyr: 12.9 ± 0.3 Gyr for the MP 
GCs, 11.8 ± 0.2 Gyr for the MR GCs. Third, the mean age for the high [a/Fe] GCs in the 
combined sample is similar to that for the low [a/Fe] GCs, while the high [a/Fe] GCs in the 
MW are, on average, older than the low [a/Fe] GCs. 

In Figure IT21 we plot the mean age and metallicity for all, MP, and MR GCs as a function 
of absolute magnitude (Mb) of gEs. The mean ages of the GCs (in the upper panels) do 
not show any systematic dependence on Mb of gEs. The mean ages of the MP GCs show 
little scatter, while those of the MR GCs show a large scatter. It is noted, however, that 
two galaxies with low luminosity (NGC 1399 and NGC 4636) have the smallest mean ages 
of MR GCs. The mean metallicities of the GCs derived from both methods (in the lower 
panels) also show little dependence on the absolute magnitudes of their host galaxies. Note 
that the results for all GCs and old GCs in gEs are not c onsistent wi t h the results derived 



from photometry of GCs for early-type galaxies in Virgo by lPeng et al.l ( 120061 ) (dashed lines) 



while those for young GCs are. Further studies are needed to explain this difference. 



5.3. [a/Fe] distribution of GCs in gEs 

Figure [131 shows the [a/Fe] distributions of the GCs in gEs. We also plot the [a/Fe] dis- 
tribution of the GCs in the combined sample for gEs and the MW GCs in panel (h). Here we 
exclude the [a/Fe] values of M87 GCs in the combined sample from the grid method because 
of abnormal value of Mg 2 Lick indices in the literature. We divide the sample into subgroups 
according to their metallicity and age as done in the previous sections. We derive the mean 
[a/Fe] of each group, which are listed in Table [TUl Notable features are as follows. First, the 
mean value of [a/Fe] for the gE GCs ( [a/Fe] = 0.14 ± 0.02 and 0.16 ± 0.01 from the grid 
method and the x 2 minimization method, respectively) is about 0.2 dex smaller than that 
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for the MW GCs, [a/Fe] = 0.36 ± 0.01. Second, the mean values of [a/Fe] for the MP and 
MR GCs in gEs are similar. Third, the mean values of [a/Fe] for the young and old GCs in 
gEs are also similar. The old GCs in the MW have much higher [a/Fe], 0.35 ±0.11 than the 
young GCs [a/Fe], —0.02 ± 0.02, although the number of the young GCs is only two. 

These [a/Fe] results of the GCs in gEs app ear to be inconsistent with those in seven 
early-type galaxies given by iPuzia et all (2005b) , wh o reported that the mean [a/Fe] for 
these gE GCs is high (0.47±0.06). IPuzia et all f l2005af ) also found that the MR GCs exhibit 
lower [a/Fe] enhancements than the MP GCs, and that [a/Fe] of GCs is independent from 
the age of GCs. The low [a/Fe] for gE GCs in this study seems to be more consistent with the 
stellar populati ons in massive ellipt ical galaxies expected from the dissipative hierarchical 



merging model ( IThomas et al.lll999l ). The discrepancy between this study and IPuzia et al. 



( I2005bl ) could be caused by the mass dif ference in their host galaxies because the our sample 
galaxies are more massive than those of IPuzia et all (j2005b| ). 



5.4. Formation of GCs in gEs seen through the metallicities and ages 



In Figures [TH and [15] we display the metallicity and age of gE GCs derived from the 
grid and x 2 minimization methods: (a) age distributions for all, MP, and MR GCs, (b) 
metallicity and age relation, and (c) metallicity distributions for all, young, and old GCs. 
The metallicities for gE GCs show an overall anti-correlation with ages, although there is a 
large scatter: the mean metallicity for GCs increases as the GCs become young. It is also 
noted that there are both MP and MR GCs in the old age, while there are only MR GCs 
in the young age. In Table [TIJ we derived the mean values of metallicity and age for the 
MR and MP GCs in the young group (age < 10 Gyr) and old group (age > 10 Gyr) in the 
combined sample of gEs and the MW. The number ratio of the young group and the old 
group is 1:3 for the MP GCs, and 1:1 for the MR GCs. In the case of the old group, the 
mean ages for the MR and MP GCs are about 13 Gyrs, and there is little difference in the 
mean ages between the MR and MP GCs. It is noted that the mean age for the MR GCs 
in the MW is about one Gyr younger than that for the MP GCs. In the case of the young 
group, the mean age for the MR GCs is about 5 Gyr, which is about 2 Gyr younger than 
that for the MP GCs. 



Puzia et all (l2005bl ) reported an existence of the age-metallicity relation based on the 



high quality spectra for 17 GCs in seven early-type galaxies. The age-metallicity relation 
shown for gE GCs as well as for MW GCs is also seen in th e studies of the GCs in the Large 
Magellanic Cloud using the high- resolution spectroscopy (jColucci et all 1201 1[). and of th e 
GCs in six nearby galaxies using the medium- resolution spectroscopy ( jSharina et al.ll2010l ). 
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Therefore, it would be very interesting to examine whether there is a global age-metallicity 
relation of GCs from dwarf galaxies to massive galaxies. 



Recently, iMuratov fc Gnedinl (120101 ) presented a semi-analytic model for the origin of 
the metallicity distribution of GCs using the galaxy merging history from the cosmological 
simulations coupled with observed scaling relations for the amount and metallicity of cold 
gas available for the star formation. They predicted that early mergers of small galaxies 
create exclusively MP GCs, while subsequent mergers with more massive galaxies create 
both MP and MR GCs. Thus they showed the age-metallicity relation, which is that MR 
GCs are several Gyr younger than MP GCs. They also expected that the mergers with more 
massive galaxies like the case of gEs would produce comparable numbers of MR and MP 
GCs simultaneously. Our results, that the MR GCs in gEs are on average ~ 3 Gyr younger 
than the MP GCs and the fra ction of young MR GCs is non-negligible, seem to be consistent 
with the model prediction by lMuratov fc Gnedinl (|2010| ). However, our result that there are 
comparable amounts of old MR GCs in gEs is not consistent with their model prediction 
that exclusively old MP GCs formed in the early stage. 



Yoon et all ( j2011bl ) suggested nonlinear color-metallicity relations between the colors of 
GCs and their metallicities. The shape of the metallicity distribution function (MDF) for 
GCs, characterized by a sharp peak with a metal-poor tail, indicates a continuous chemi- 
cal enrichment with a short timescale less than 1 Gyr. They further suggested a possible 
systematic age difference among GC systems, in that the GC systems in more luminous 
galaxies are older. Their result related to the MDF shape is not consistent with our result 
that the old GCs show a bimodal metallicity distribution. In addition, their result about the 
systematic age difference among GC systems is not consistent with our result in the sense 
that the mean ages of GCs in gEs do not show any systematic dependence on the luminosity 
of gEs as shown in Figure [T2"l 

Preyiously, several formation models of the GC systems in gEs have been suggested 



flPeebles 111969 : lAshman fc Zepjll992l : Harris et al.lll99a Forbes et al.lll997t ICote et al.lll998f ) 



(Rhode & ZeDf 


2001 


Lee 


2003; 


Ric 


itler et al. 


2004 


West et al. 


2004; 


Brodie & Strader 


2006; 


Hwang et al. 


2008; 


Lee et al.l 


2010a 


) . Considering the observational results on the kinematics 



and photometry of GCs in gEs and predictions of several models available in the literature, 
Lee et al.l (l2010al ) presented a mixture (bibimbap) scenario describing the formation of GCs 
in gEs. However, they did not use the information on chemical properties and ages of GCs in 
gEs when they proposed their scenario. To be short, the scenario is as follows. (1) MP GCs 
are formed mostly in low-mass dwarf galaxies very early, and preferentially in dwarf galaxies 
located in the high density environment like galaxy clusters. These are the first generation 
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of GCs in the universe. MP GCs should be also formed in massive galaxies as well, but the 
number of these massive galaxies is much smaller than that of the dwarf galaxies. (2) MR 
GCs are formed together with stars in massive galaxies or dissipational merging galaxies 
later than MP GCs, but not much later than MP GCs. The chemical enrichment of the 
galaxies is rapid after the formation of MP GCs and the difference in the formation epoch 
of MP GCs and MR GCs should be small. (3) Massive galaxies grow becoming gEs via 
dissipationless or dissipational merging of galaxies of various types and via accretion of 
many dwarf galaxies. New MR GCs will be formed during the dissipational merging, but 
the fraction of dissipational merging at this stage should be minor. A significant fraction of 
MP GCs in gEs we see today are from dissipationless merging or accretion. 

If we consider our results on ages and chemical properties in the view of thi s mixture 



scenari o, we can conclude the followings. First, the point in the mixture scenario of iLee et al. 



( 12010al ) that MP GCs in massive galaxies are formed early is consistent with our finding that 
most MP GCs in gEs are older than 10 Gy r. Second, the sm all difference in the primary 
formation epochs of MP GCs and MR GCs in lLee et al.l (l2010al ) is consistent with our results 
that there is little difference in the mean ages of the MP and MR sub-populations of old 
GCs (> 10 Gyr). Note that there is ~ 1 Gyr difference in the case of the MW GCs that 
have more reliable age estimates (12.9 ± 0.3 Gyr and 11.8 ± 0.2 Gyr for the MP and MR 
GCs, respectively). 

Important findings in this study are that gEs have a significant number of young GCs 
and that they are mostly MR GCs. We also found that there are a smaller number of 
MP GCs with young ages (< 10 Gyr). The MR young GCs are probably formed during 
the dissipational merging of late-type galaxies with their host elliptical galaxies, while the 
MP young GCs are formed in dwarf galaxies and are accreted to their host galaxies. It is 
expected that the fraction of new MR GCs from the dissipational merging events depends 
on the environment such as t he gas content of the merging compani ons and the local density 
where the merging happens (IPark fc Choill2009i ; iHwang et al.ll2010l ). Therefore, the fraction 
of young GCs would not be the same for all galaxies. In the case of NGC 4636, a considerable 
fraction of young GCs with ages less than 5 Gyr might be the result of recent merging with 
gaseous late-type galaxies. It would be worth investigating the environmental dependence of 
the GC systems in gEs once the number of gEs with spectroscopically measured metallicity 
and age for their GCs increases in the future. 
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6. Summary 

We presented the measurements of the metallicity, age, and [a/Fe] of GCs in NGC 4636 
derived from the Subaru spectroscopic data. For comparison we have also re-derived the 
metallicity, age, and [a/Fe] of GCs for six other gEs (M87, M49, M60, NGC 5128, NGC 
1399, and NGC 1407) from the line indices in the literature using the same methods as used 
for NGC 4636 GCs. We made a combined sample of GCs in gEs by combining all the GC 
data in gEs including NGC 4636, and used it to investigate various properties of GCs in 
gEs. Final results for the combined sample as well as individual galaxies are summarized in 
Tables 7-11. Our primary results are summarized below. 

1. We measured the metallicities of 59 GCs in NGC 4636 with the BH method, and the 
metallicity, age, and [a/Fe] values of 33 GCs with the grid method and x 2 minimization 
method. 

2. The metallicities of the NGC 4636 GCs show marginally a bimodal distribution with 
two peaks at [Fe/H] = -1.23(cr = 0.32) and -0.35(o = 0.19) from the BH method, 
-1.14(a = 0.22) and -0.41(o = 0.21) from the grid method, and -1.22(o = 0.25) and 
— 0.44(o = 0.16) from the % 2 minimization method. The mean values of metallicity 
are [Fe/H] = —0.90 ± 0.12 from the BH method, —0.70 ± 0.13 from the grid method, 
and —0.74 ±0.11 from the x 2 minimization method. 

3. The ages of the GCs in NGC 4636 show a large range with a non-negligible number 
of young GCs. The mean age of the GCs is derived to be 7.9 ± 0.8 Gyr from the grid 
method and 8.7 ± 2.0 Gyr from the x 2 minimization method. The number ratio for 
young (< 5 Gyr), intermediate (5 — 10 Gyr), and old (> 10 Gyr) GCs is 27:39:34 from 
the grid method and 27:30:43 from the x 2 minimization method. 

4. The [a/Fe] distribution of the NGC 4636 GCs is broad with the mean of [a/Fe] = 
0.20 ±0.11 from the grid method and of [a/Fe] = 0.07±0.05 from the x 2 minimization 
method, smaller than the value for the MW GCs. 

5. The metallicity, age, and alpha-element of the GCs in NGC 4636 do not show any 
significant radial variation. NGC 4636 GCs show a hint of the age-metallicity relation. 

6. The GCs in the combined sample of gEs show a bimodal metallicity distribution: 
peaks at [Fe/H] = — 1.21(cr = 0.32), and -0.42(o = 0.27) from the grid method, and 
[Fe/H] = -1.25(<T = 0.32), and -0.42(o = 0.25) from the % 2 minimization method. 
The MP GCs are on average ~ 3 Gyr older than the MR GCs. They show a large 
range in age from 2 to 15 Gyr, including a significant number of young GCs, while all 
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MW GCs are older than 7 Gyr. The MR GCs show a broader age distribution than 
the MP GCs, showing that young GCs are mostly metal-rich. The mean [a/Fe] value 
of the combined gE GCs is smaller than that of the MW GCs. 

7. The number ratio of the young group and the old group is 1:3 for the MP GCs, and 
1:1 for the MR GCs. In the case of the old group, the mean ages for the MR and MP 
GCs are about 13 Gyrs, and there is little difference in the mean ages between the MR 
and MP GCs. It is noted that the mean age for the MR GCs in the MW is about one 
Gyr younger than that for the MP GCs. In the case of the young group, the mean age 
for the MR GCs is about 5 Gyr, which is about 2 Gyr younger than that for the MP 
GCs. 

8. The results on metallicity, age, and [a/Fe] of GCs in gEs as well as t hose on ph o tometr y 



and kinematics of GCs are consistent with the mixture scenario in iLee et all ( l2010al ). 
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Table 1. Zero Point Offsets for Lick Index calibration 



index 


unit 


onset 


rms 


CN1 


mas' 


-0.021 


0, 


004 


CN2 


mas' 


-0.007 


0. 


.018 


Ca4227 


A 


0.021 


0. 


.367 


G4300 


A 


0.414 


0. 


.249 


Fe4383 


A 


0.515 


0. 


.259 


Ca4455 


A 


0.408 


0. 


.291 


Fe4531 


A 


0.842 


0. 


.780 


C24668 


A 


-1.880 


1. 


036 


H/3 


A 


0.184 


0. 


.140 


Fe5015 


A 


0.616 


0. 


.267 


Mgl 


mag 


-0.017 


0. 


016 


Mg2 


mag 


-0.015 


0. 


.018 


Mgb 


A 


0.165 


0. 


.187 


Fe5270 


A 


0.443 


0. 


.306 


Fe5335 


A 


-0.103 


0. 


.162 


Fe5406 


A 


0.271 


0. 


.147 


Fe5709 


A 


0.148 


0. 


.484 


Fe5782 


A 


0.036 


0. 


.156 


NaD 


A 


0.494 


0. 


.578 




A 


-0.367 


0. 


.458 


H 7 A 


A 


-0.761 


0. 


.645 


m F 


A 


0.129 


0. 


.242 


H 7 F 


A 


0.083 


0. 


269 



Index(Lick)=Index(Sub&ru)+constant. 
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Table 2. Lick Line Indices (set 1) 



ID 11 


CN1 


CN2 


Ca4227 


G4300 


Fe4383 


Ca4455 


Fe4531 


C24668 


H/3 


Fe5015 


Mgl 


Mg2 




(mag) 


(mag) 


(A 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(mag) 


(mag) 


157 


0.219 


0.242 


0.747 


4.612 


6.693 


-0.054 


4.808 


-2.320 


1.865 


4.851 


0.093 


0.219 


182 


-0.070 


-0.030 


1.912 


0.824 


-1.213 


0.581 


0.719 


-3.783 


1.935 


-0.530 


0.001 


0.073 


244 


0.290 


0.322 


-3.348 


6.267 


5.615 


3.151 


1.416 


-0.785 


2.230 


3.279 


0.056 


0.233 



a From iPark et all j201& |2012a|) 



Note. — This tabic is available in its entirety in amachinc-rcadable form in the online journal. A portion is shown here for guidance 
regarding its form and content. 



Table 3. Lick Line Indices (set 2) 



ID" 


Mgb 


Fe5270 


Fe5335 


Fe5406 


Fe5709 


Fe5782 


NaD 


HS A 


H7A 


H<5 F 


H7j7 




(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


157 


3.837 


2.943 


2.952 


2.993 


1.325 


0.556 


1.426 


-3.797 


-5.303 


-0.301 


-0.951 


182 


0.322 


1.199 


2.723 


-0.436 


0.916 


0.360 


5.777 


-0.220 


3.506 


-0.837 


2.851 


244 


5.220 


0.723 


2.289 


1.335 


1.666 


-0.001 


1.963 


-3.485 


-17.252 


-1.465 


-7.051 



a From Park et all <201Cl 12012^ 

Note. — This table is available in its entirety in amachinc-rcadable form in the online journal. A portion is shown here 
for guidance regarding its form and content. 



Table 4. Lick Index Errors (set 1) 



ID a 


CN1 


CN2 


Ca4227 


G4300 


Fe4383 


Ca4455 


Fo4531 


C24668 


H/3 


Fe5015 


Mgl 


Mg2 




(mag) 


(mag) 


(A 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(mag) 


(mag) 


157 


0.024 


0.031 


0.502 


0.931 


1.217 


0.553 


0.898 


1.341 


0.447 


1.013 


0.014 


0.016 


182 


0.019 


0.025 


0.404 


0.790 


1.188 


0.502 


0.916 


1.367 


0.439 


1.019 


0.014 


0.016 


244 


0.028 


0.036 


0.614 


1.168 


1.513 


0.633 


1.082 


1.574 


0.502 


1.143 


0.016 


0.012 



a From IPark et all j201Cl |2012a|) 



Note. — This table is available in its entirety in amachinc-rcadable form in the online journal. A portion is shown here for guidance 
regarding its form and content. 



Table 5. Lick Index Errors (set 2) 



ID a 


Mgb 


Fe5270 


Fe5335 


Fe5406 


Fe5709 


Fe5782 


NaD 


HS A 


H 7 A 


HS F 


H 7F 




(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


157 


0.457 


0.461 


0.581 


0.394 


0.275 


0.266 


0.424 


1.013 


0.756 


0.647 


0.438 


182 


0.435 


0.454 


0.567 


0.409 


0.275 


0.278 


0.423 


0.848 


0.655 


0.562 


0.397 


244 


0.516 


0.525 


0.678 


0.444 


0.308 


0.289 


0.437 


1.148 


0.918 


0.751 


0.508 



a From Park et al. ( 201(1 12012^ 



Note. — This table is available in its entirety in amachinc-rcadable form in the online journal. A portion is shown 
here for guidance regarding its form and content. 
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Table 6. Metallicity, Age, and [a/Fe] of NGC 4636 GCs 



ID a [Z/H] lM . 3 / e] ' 


[Z/H] M 9 2 




[a/Fe] gr i d 


[Z/H] X 2 


Age x 2 


[a/Fe] % 2 


[Fe/H] SH 


(dex) 


(dex) 


(Gyr) 


(dex) 


(dex) 


(Gyr) 


(dex) 


(dex) 



157 0.27 ±0.28 0.11 ±0.05 4.7 ±3.9 -0.08 ±0.10 0.11 ±0.09 7.0 ± 2.6 -0.20 ±0.12 -0.13 ±0.16 
182 -1.70 ±0.39 -1.04 ±0.10 14.0 ±1.5 -0.30 ± 0.00 -1.46 ±0.23 12.0 ± 3.6 -0.30 ± 0.25 -1.79 ±0.16 
244 -0.21 ±0.22 -0.12 ±0.05 5.4 ± 4.0 0.49 ± 0.02 0.23 ± 0.09 2.2 ± 0.8 0.26 ± 0.21 -0.32 ± 0.48 



a From lPark et ail j20ld l2012ah 

Note. — This table is available in its entirety in amachine-readable form in the online journal. A portion is shown here for 
guidance regarding its form and content. 



Table 7. KMM Test Results for the Metallicities of GCs in gEs ; 



galaxy 


N 


[Fe/H] mp ± cr 
(dex) 


N 


[Fe/H] M fl ± <x 
(dex) 


confidence level 
(%) 






grid method 






NGC 4636 


10 


-1.14 ±0.22 


15 


-0.41 ± 0.21 


71.4 


NGC 5128 


12 


-1.14 ±0.32 


11 


-0.34 ±0.12 


97.7 


M60 


5 


-1.27 ±0.62 


33 


-0.56 ± 0.36 


71.9 


NGC 1407 


7 


-1.14 ±0.23 


12 


-0.38 ± 0.18 


76.2 


NGC 1399 


3 


-1.63 ± 0.06 


7 


-0.25 ± 0.51 


94.6 


M87 


66 


-1.24 ±0.32 


63 


-0.45 ± 0.27 


97.5 


M49 


28 


-0.82 ± 0.65 


15 


0.37 ±0.10 


99.9 


gEs 


155 


-1.21 ± 0.32 


119 


-0.42 ± 0.27 


99.9 


gEs (< 10 Gyr) 


22 


-1.41 ± 0.14 


94 


-0.57 ±0.37 


98.4 


gEs (> 10 Gyr) 


111 


-1.22 ± 0.34 


77 


-0.42 ± 0.25 


99.7 






y 2 minimization method 




NGC 4636 


14 


-1.22 ±0.25 


14 


-0.44 ±0.16 


95.2 


NGC 5128 


13 


-1.40 ±0.47 


19 


-0.43 ± 0.15 


98.2 


M60 


20 


-1.09 ± 0.31 


15 


-0.37 ±0.11 


98.8 


NGC 1407 


8 


-1.05 ± 0.33 


11 


-0.39 ± 0.11 


92.8 


NGC 1399 


4 


-1.68 ± 0.24 


6 


-0.37 ±0.21 


94.5 


M87 


40 


-1.41 ± 0.14 


90 


-0.67 ±0.36 


99.9 


M49 


26 


-1.15 ± 0.71 


20 


0.19 ±0.25 


99.2 


gEs 


167 


-1.25 ± 0.32 


148 


-0.42 ± 0.25 


99.9 


gEs (< 10 Gyr) 


63 


-1.11 ± 0.43 


64 


-0.32 ± 0.27 


98.5 


gEs (> 10 Gyr) 


114 


-1.25 ± 0.29 


74 


-0.44 ±0.18 


99.9 


MW b 


111 


-1.52 ± 0.39 


41 


-0.52 ± 0.23 


99.8 



a Peak values, widths (a), and confidence levels are derived from the KMM test. 
b It is from lHarrisI l|l996h . 



Table 8. Mean Metallicities of GCs in gEs a 







total 




age < 10 Gyr 




age > 10 Gyr 




[a/Fe] < 0.2 




[a/Fe] > 0.2 


galaxy 


N 


<[Fe/H]> ± < a > 


N 


<[Fe/H]> ± < cr > 


N 


<[Fe/H]> ± < cr > 


N 


<[Fe/H]> ± < cr > 


N 


<[Fe/H]> ± < cr > 






(dex) 




(dex) 




(dex) 




(dex) 




(dex) 












grid method 










IN 40o0 


oo 


n 7Pi _i_ n i q 
— U. / U ± U.lo 


9 9 

zz 


n a o -l- pi 1/t 


11 


-1.20 ±0.25 


"1 7 


— U.OU it U.lo 


10 


O _|_ n 91 
— U.oO it U..Z1 


IN \J\^ DlZo 


f-: o 


i in -l- n os 

— 1 . 1U It U.UO 


zo 


— U.ol it U.±l 


45 


-1.27±0.11 


11 


1 07 4-0 11 
— 1 .u / It U. 1 1 


9/1 
Z 1 


i 1 a A- n 1 ^ 

— ±.144 It U. IO 


1V10U 


oo 


o fi7 -i- n os 

— U.O / It U.UO 


i n 

1U 


n Q'l. 4- o 90 
— u.yo it u.iU 


28 


-0.59 ±0.09 


9-1 
Z'l 


70 4- O OQ 

— u. (u it u.uy 


"1 A 
11 


O ^£ 4- O 1 S 
— U.Oo It U. lo 


WPP i MY7 




— U.Oo It U. lo 


i 
i 


— U.4o It u.zu 


16 


-0.74 ±0.15 


1 

1 


9 1Q 4- O 1 fi 

— z.oy it u. io 


iy 


n fip; 4-01^ 
— u.oo It u. io 


MPP 1 QQQ 

in ioyy 


"1 


pi _|_ pi 97 


7 


n 'zi 4- o 47 

— U.OO It U.4t ( 


3 


-1.60 ±0.43 


7 


4- O AQ 

— u.yo it u.4y 


Q 
O 


O 39 4- O f}9 
U.OZ It u.oz 


IVlo ( 


1 A (\ 
1 10 


O fiQ 4- O C\R 

— U.oo it U.UO 


OO 


n ^7 -U n HQ 

— U.O / it U.Uo 


90 


-1.00 ±0.07 










M49 


43 


-0.42 ±0.14 


20 


-0.10 ± 0.34 


23 


-0.68 ±0.20 


30 


-0.48 ±0.16 


13 


-0.22 ±0.30 


gEs 


358 


-0.80 ±0.04 


142 


-0.55 ± 0.05 


216 


-0.96 ±0.05 


123 


-0.78 ±0.08 


89 


-0.76 ±0.08 










x 2 


minimization method 










NGC 4636 


33 


-0.74 ±0.11 


19 


-0.59 ±0.16 


14 


-0.94 ±0.14 


20 


-0.92 ±0.13 


13 


-0.40 ±0.11 


NGC 5128 


62 


-1.15 ±0.08 


21 


-1.11 ± 0.14 


41 


-1.17 ±0.09 


16 


-1.20 ±0.10 


16 


-1.01 ±0.14 


M60 


38 


-0.76 ±0.09 


11 


-0.58 ± 0.26 


27 


-0.82 ±0.11 


16 


-0.75 ±0.15 


22 


-0.76 ±0.13 


NGC 1407 


20 


-0.62 ±0.13 


6 


-0.52 ± 0.16 


14 


-0.69 ±0.17 


6 


-0.91 ±0.19 


14 


-0.44 ±0.12 


NGC 1399 


10 


-0.83 ±0.40 


7 


-0.51 ± 0.36 


3 


-1.16 ±0.79 


5 


-1.48 ±0.78 


5 


-0.56 ±0.20 


M87 


150 


-0.90 ±0.05 


52 


-0.63 ± 0.08 


98 


-1.08 ±0.06 


77 


-1.02 ±0.07 


73 


-0.77 ±0.08 


M49 


46 


-0.58 ±0.16 


24 


-0.44 ±0.18 


22 


-0.77 ±0.25 


23 


-0.69 ±0.19 


23 


-0.38 ±0.29 


gEs 


359 


-0.86 ±0.04 


140 


-0.65 ± 0.06 


219 


-1.00 ±0.04 


193 


-1.00 ±0.05 


166 


-0.70 ±0.05 


MW 


152 


-1.28 ±0.05 


3 


-0.61 ± 0.27 


61 


-1.44 ±0.08 


8 


-1.33 ±0.30 


49 


-1.40 ±0.10 



to 

CO 



a The mean values are determined with the biweight location of lBeers et al. The uncertainties, < cr >, indicate the 68% confidence levels obtained with a 

bootstrap method. 



Table 9. Mean Ages of GCs in gEs : 







total 




[Fe/H] < -0.9 




[Fe/H] > -0.9 




[a/Fe] < 0.2 




[a/Fe] > 0.2 


galaxy 


N 


<age> ± < a > 


N 


<age> ± < <t > 


N 


<age> ± < a > 


N 


<age> ± < ct > 


N 


<age> ± < a > 
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grid method 
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M49 


43 


9.5 ±2.1 


13 


12.1 ± 0.2 


30 


7.7 ± 1.2 


30 


8.0 ± 1.3 


13 


12.2 ± 0.6 


gEs 


358 


12.3 ± 0.1 


157 


12.2 ± 0.1 


201 


9.2 ±0.7 


123 


12.4 ± 1.5 


89 


12.2 ± 0.2 










x 2 


minimization method 










NGC 4636 


33 


8.7 ±2.0 


14 


11.9 ± 1.5 


19 


7.3 ±1.2 


20 


9.0 ± 1.0 


13 


8.2 ±3.4 


NGC 5128 


(52 


12.1 ± 0.4 


10 


12.3 ± 0.3 


22 


11.9 ± 1.3 


16 


12.4 ±0.5 


16 


11.5 ± 1.0 


M60 


38 


12.9 ± 0.3 


16 


12.6 ± 0.4 


22 


13.0 ±0.6 


16 


12.7 ± 1.3 


22 


12.8 ± 0.2 


NGC 1407 


20 


12.8 ±0.5 


6 


12.9 ±0.3 


14 


12.6 ± 1.4 


6 


10.0 ±2.1 


14 


13.0 ±0.2 


NGC 1399 


10 


7.5 ± 1.8 


•1 


11.5 ± 2.0 


6 


5.1 ±1.7 


5 


9.5 ±3.3 


5 


4.8 ±2.6 


M87 


150 


12.3 ±0.1 


75 


12.2 ± 0.1 


75 


11.7 ±2.1 


77 


12.3 ± 0.1 


73 


12.4 ±0.1 


M49 


46 


9.5 ±1.5 


18 


11.7 ± 1.0 


28 


8.6 ±1.2 


23 


9.0 ± 1.0 


23 


11.8 ± 1.5 


gEs 


359 


12.2 ± 0.1 


173 


12.2 ± 0.1 


186 


9.7±0.7 


193 


12.2 ± 0.1 


166 


12.3 ± 0.1 


MW 


64 


12.2 ± 0.1 


-18 


12.9 ± 0.3 


16 


11.8 ± 0.2 


6 


10.7 ± 1.0 


38 


13.0 ±0.2 



CO 

o 



a The mean values are determined with the biweight location of Be ers et al. I (TjjjgiJ). The uncertainties, < a >, indicate the 68% confidence levels 
obtained with a bootstrap method. 



Table 10. Mean [a/Fe]s of GCs in gEs : 



^a,! ax v 


N 


total 

<[a/Fe] > ± < a > 
(dex) 


N 


[Fe/H] < -0.9 
<[a/Fe] > ± < a > 
(dex) 


N 


[Fe/H] > -0.9 
<[a/Fe] > ± < a > 
(dex) 


N 


age < 10 Gyr 
<[a/Fe] > ± < a > 
(dex) 


N 


age > 10 Gyr 
<[ct/Fe] > ± < a > 
(dex) 












grid method 










NGC 4636 


33 


0.20 ± 0.11 


13 


0.48 ± 0.15 


20 


0.17 ±0.06 


22 


0.22 ± 0.14 


11 


0.17 ± 0.20 


NGC 5128 


68 


0.13 ± 0.03 


42 


0.14 ± 0.04 


26 


0.12 ±0.05 


23 


0.14 ± 0.05 


15 


0.12 ± 0.04 


M60 


38 


0.06 ± 0.04 


12 


0.01 ± 0.22 


26 


0.08 ±0.05 


10 


-0.08 ± 0.15 


28 


0.09 ± 0.05 


NGC 1407 


20 


0.42 ± 0.05 


8 


0.44 ± 0.05 


12 


0.40 ± 0.06 


1 


0.36 ± 0.06 


16 


0.46 ± 0.05 


NGC 1399 


10 


0.07 ± 0.09 


3 


-0.30 ± 0.19 


7 


0.21 ±0.08 


7 


0.09 ± 0.10 


3 


0.09 ± 0.39 


M87 






















M49 


43 


0.09 ±0.04 


13 


-0.04 ±0.19 


30 


0.13 ±0.05 


20 


0.02 ± 0.05 


23 


0.15 ±0.07 


gEs 


212 


0.14 ±0.02 


91 


0.14 ±0.03 


121 


0.15 ±0.02 


86 


0.12 ± 0.03 


126 


0.16 ±0.02 










x 2 


minimization method 










NGC 4636 


33 


0.07 ±0.05 


14 


-0.10 ±0.07 


19 


0.21 ±0.07 


19 


0.05 ± 0.06 


14 


0.10 ±0.11 


NGC 5128 


62 


0.06 ±0.03 


10 


0.02 ±0.04 


22 


0.12 ±0.03 


21 


0.07 ±0.06 


41 


0.05 ±0.03 


M60 


38 


0.26 ±0.05 


16 


0.36 ±0.15 


22 


0.26 ±0.03 


11 


0.24 ±0.08 


27 


0.28 ±0.05 


NGC 1407 


20 


0.31 ±0.03 


6 


0.21 ±0.13 


11 


0.33 ±0.04 


6 


0.25 ± 0.07 


14 


0.33 ±0.04 


NGC 1399 


10 


0.15 ±0.09 


1 


-0.25 ± 0.15 


6 


0.32 ±0.15 


7 


0.18 ± 0.19 


3 


0.17 ±0.18 


M87 


150 


0.16 ±0.02 


75 


0.05 ± 0.04 


75 


0.24 ±0.02 


52 


0.21 ± 0.04 


98 


0.13 ±0.03 


M49 


16 


0.19 ±0.04 


18 


0.11 ±0.14 


28 


0.24 ±0.04 


24 


0.17 ±0.06 


22 


0.23 ±0.11 


gEs 


359 


0.16 ±0.01 


173 


0.07 ±0.02 


186 


0.23 ±0.01 


140 


0.16 ± 0.02 


219 


0.15 ±0.02 


MW 


57 


0.36 ±0.01 


43 


0.38 ±0.01 


14 


0.29 ±0.04 


2 


-0.02 ±0.01 


12 


0.37 ±0.01 



a The mean values are determined with the biweight location of lBeers et al ] l|l990h . The uncertainties, < a >, indicate the 68% confidence levels obtained with a 
bootstrap method. 
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Table 11. Mean Age and Metallicity of GCs in gEs and the MW' 



young GC (age < 10 Gyr) 
N <[Fe/H]> <age> 
(dex) (Gyr) 



old GC (age > 10 Gyr) 
N <[Fe/H]> <age> 
(dex) (Gyr) 



gEs (grid) 

MP GC ([Fe/H] < -0.9) 42 -1.26 ±0.04 7.0 ±0.3 115 -1.41 ±0.06 12.5 ± 0.4 

MR GC ([Fe/H] > -0.9) 100 -0.27 ±0.05 4.8 ±0.3 101 -0.45 ± 0.03 13.0 ± 0.3 

gEs (y 2 minimization) 

MP GC ([Fe/H] < -0.9) 50 -1.33 ±0.05 7.2 ±0.3 123 -1.39 ±0.06 12.5 ± 0.5 

MR GC ([Fe/H] > -0.9) 90 -0.32 ±0.04 5.1 ±0.3 96 -0.48 ± 0.03 13.0 ± 0.3 

MW 

MP GC ([Fe/H] < -0.9) ... ... 48 -1.64 ±0.05 12.9 ± 0.3 

MR GC ([Fe/H] > -0.9) 3 -0.61 ± 0.26 7.0 ± 0.8 13 -0.60 ± 0.06 11.8 ±0.3 



a The mean values are determined with the biweight location of iBeers et all (ll990f ). The 
uncertainties indicate the 68% confidence levels obtained with a bootstrap method. 
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Fig. 1.- Example spectra of (a) a metal-rich (MR) GC in NGC 4636 (ID: 154) with 
T\ = 19.82 and [Fe/H] = -0.39 dex; (b) a metal-poor (MP) GC in NGC 4636 (ID: 9995) 
with 7\ = 19.22 and [Fe/H] = -1.49 dex; (c) NGC 4636 nucleus; (d) NGC 6624, a MR MW 
GC with [Fe/H] = -0.44 dex; and (e) M13, a MP MW GC with [Fe/H] = -1.54 dex. All 
spectra are plotted in the rest frame and are smoothed with the Lick resolution. 
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Fig. 2. — Lick line index diagrams, (a) H/3 line versus [MgFe]' fo r NGC 4636 GCs me asured 
in this study (circles). The crosses represents NGC 5128 GCs ( iWoodley et al.ll2010l ). The 
grids represent the SSP models for various values of [Z/H] (-2.25, -1. 35, -0.33, 0, 0.35, and 
0.67) and ages (0.4, 0.6, 0.8, 1, 2, 3, 5, 8, 10, and 15 Gyr) given bv iThomas et all tooA 



20041 ). (b) <Fe> versus Mg 2 . The grids are for a model age of 8 Gyr. 



-3 -2-10 1 -3 -2-10 1 -3 -2-10 1 
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-3 -2 -1 1 5 10 15 -0.5 0.5 
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Fig. 3. — (a) Comparison of metallicities for NGC 4636 GCs measured with the BH method 
and with the grid method. Here the grid metallcities are the values derived from the B./3 ver- 
sus [MgFe]' grid, (b) Comparison of [Fe/H] for NGC 4636 GCs measured with the BH 
method and the x 2 minimization method (c) Comparison of [Z/H] values obtained from 
Mg 2 grids and [MgFe]' grids. Comparison of parameters for NGC 4636 GCs measured with 
the grid method and the x 2 minimization method: (d) [Z/H], (e) age, and (f) [a/Fe]. Here 
the grid metallcities are the values derived from the H/3 versus [MgFe]' grid. The solid 
and dotted lines represent the linear least-squares fit with 2 a clipping and the one-to-one 
relation, respectively. 
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Fig. 4. — Metallicity, age, and [a/Fe] distribution for NGC 4636 GCs. (a) Distribution of 
[Fe/H] from the BH method, (b) [Z/H] distributions obtained from H/3 versus [MgFe]' grid, 
(c) Age distribution, (d) [a/Fe] distribution. The solid and dot-dashed histograms represent 
the distributions derived with grid and \ 2 method, respectively. The error bars in the upper 
right of each panel indicate the mean error of each value. 
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Fig. 5.— Metallicity distribution of NGC 4636 GCs. (a) [Fe/H] distribution from t he BH 
meth od. The hatched histogram represents the [Fe/H] distribution of the MW GCs ([Harris 
19961 ). The histograms in (b) and (d) represent [Fe/H] distributions obtained from the grid 
and x 2 minimizati on method, respect ively, (c) Photometric metellicity d istribution o btaine d 
from (C — Ti)o by IPark et al.l (j2012al ) and the double linear equation by lLee et al. I (]2008al ). 
The curved solid lines represent double Gaussian fits. 
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Fig. 6. — Radial variation of metallicity, age, and [a/Fe]: (a) [Fe/H] from the BH method, 
(b) [Z/H], (c) Age, and (d) [a/Fe] from the grid method, and (e) [Z/H], (f) Age, and (g) 
[a/Fe] from the x 2 method. The solid and dotted lines represent the mean values and their 
dispersions derived using a moving radial bin (with 2.5 arcmin width and 0.5 arcmin step). 
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Fig. 7. — Relations between metallicity and age from the grid (a) and the x 2 minimization 
method (b). The circles and cross es represent the NGC 4636 GCs in this study and the NGC 
5128 GCs derived from indices in IWoodlev et al. respectively. The squares represent 

the MW GCs by lHarrisI (119961 ) and ICarretta et al.l (120101 ) . The large boxes represent the 
mean errors for NGC 4636 GCs. 



-40- 



8 
6 
' 4 

2 


10 

8 

6 

4 

2 






& 2 



10 
8 
6 
4 
2 




n — i — i — | — i — i — i — — | — i — r 

(a) NGC 46p6 f."; 



(c) M60 



I I I I I I; I 



— i — i — | — i — i — r 
age^ lOGyr- 
age> lOGyr 



1 1, I I I I 



4 -(e) NGC 1^99 



M l 



"(g) M49 




I M, I I I 



L.J 



M 



J I I I I I I I I I L 



-2 -1 
[Fe/H] 



10 - 




6 




25 
20 
15 
10 
5 


40 
20 




(b) NGC 5128 i 



I !r- 
U-i 



l_J ■ ■" i "■ " ' 



M ill 



(d) NGC 1407 



I I I I 



n 



(f) M87 



(h) 7 gEs 
grid 



M rr' 



l.j 



j i i i i i_ 



J I I I L_ 



-2 -1 
[Fe/H] 



Fig. 8. — Metallicity distribution of GCs in gEs from the grid method. The solid, dot- 
dashed, and dotted line histograms for gEs show the [Fe/H] distributions of all GCs, old 
GCs (> 10 Gyr), and young GCs (< 10 Gyr), respectively. The vertical dotted and dot- 
dashed lines represent the mean metallicities for young and old GCs, respectively. In panel 
(h) the histograms represent the metallicity distribution for the combined sample of GCs in 
all seven gEs. 
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Fig. 9. — Metallicity distribution of GCs in gEs from the x 2 minimization method. Symbols 
are same as Figure [HJ 
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Fig. 10. — Age distribution of the MP ([Fe/H] < —0.9) GCs (shaded histograms) and the 
MR ([Fe/H] > —0.9) GCs (hashed histograms) in gEs derived from the grid method. The 
solid line histograms represent the sum of the MP and MR groups. The mean ages of the MP 
and MR groups are shown by the dotted and dot-dashed lines, respectively. The histograms 
on panel (h) show the age distributions for the combined sample of GCs in seven gEs. 
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Fig. 11.— Age distribution of GCs in gEs derived from the \ 2 minimization method. 
Symbols are same as Figure [TDJ 
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Fig. 12. — Mean ages and metallicities of GCs as a function of absolute magnitude (Mb) 
of gEs. (a,d) All GCs, (b) MP GCs ([Fe/H] < -0.9), (c) MR GCs ([Fe/H] > -0.9), (e) 
young GCs ( age < 10 Gyr), and (f) old GCs (age > 10 Gyr). M B 's of gEs are obtained from 
HyperLeda jPaturel et aljboO.sl ). The filled and open symbols are from the grid method 
and the x 2 minimization method, respectively. The dashed lines represent a r elation derived 
from photometry of GCs in early-type galaxies in Virgo by iPeng et al.l (120061 ). 
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Fig. 13.— [a/Fe] distribution of GCs in gEs from the grid method (solid) and the x 2 
minimization method (dot-dashed). The vertical lines represent the mean values of the 
[a/Fe] of GCs in each galaxy. Panel (h) shows the [a/Fe] distribution for the combined 
sample of GCs in six gEs (grid method) excluded M87 GCs and in seven gEs (x 2 fit). The 
dotted histogram and line represent the MW GCs. 
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Fig. 14. — Age-metallicity relation of GCs in gEs from the grid method, (a) The solid, 
shaded, and hashed histograms represent the age distributions of all, MP, and MR GCs, 
respectively, (c) The dotted and dot-dashed histograms show the [Fe/H] distributions for 
the combined sample of GCs with ages smaller and larger than 10 Gyr, respectively, in seven 
gEs. The solid histogram represents the metallicity distribution of all GCs. (b) The squares, 
circles, triangles, pentagons, diamonds, plus, and stars represent the GCs of M60, NGC 
4636, NGC 5128, NGC 1407, NGC 1399, M87, and M49, respectively. The vertical dotted 
line is an age criterion, 10 Gyr, to divide young and old GCs and the horizontal dashed line 
is a metallicity criterion, [Fe/H] = —0.9 dex, to divide MP and MR GCs. 
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Fig. 15. — Age-metallicity relation of GCs in gEs from the x 2 minimization method. Symbols 
are same as Figure [TH 



